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Abstract 

Heavy quark production in five jet events at e+e~ colliders is studied 
at tree level using helicity amplitudes. Total production rates for 263j and 
Ahj are given and compared with the corresponding results for massless 
quarks. The process e+e~ — >■ ggggj which is the dominant contribution to 
4j7 production is briefly discussed. 
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Introduction 



The great number of hadronic decays of the Z observed at LEP provides the oppor- 
tunity to test our understanding of strong interactions in unprecedented detail. Large 
samples of multi-jet events have been accumulated and analyzed [0, 0. In particular 
the ALEPH Collaboration and the OPAL Collaboration have studied the jet-fraction 
distributions as a function of i/cut for up to five jets. 

Recent advances in 6-tagging techniques based on the introduction of vertex detec- 
tors and on a refinement of the selection procedures, with their large efficiencies and 
the resulting high purities, have paved the way to the study of heavy-quark production 
in association with light-quark and gluon jets. The flavour independence of the strong 
coupling constant has been verified by several groups p using both jet-rates and shape 
variables. 

It has been lately pointed out |^ that the effects of the 6-quark mass are sub- 
stantial and increase with the number of jets. 

In this paper we plan to study heavy-quark production in five-jet events at e~^e~ 
colliders. Two different mechanisms contribute to five-jet production for light quarks, 
the 'point-like' annihilation of e~^e~ to Z^, 7 |Q, which dominates at LEP I energy, and 
W^W~ production [0, followed by the emission of a gluon from the decay products of 
the VT's, which dominates above the WW treshold. This latter contribution is severely 
suppressed for 6-quarks by the smallness of the Vhc element of the CKM matrix and 
by the large mass of the top and will not be considered here. We have studied the 
annihilation mechanism for five-jet production at tree level taking full account of 7, Z 
interference and of quark masses. Therefore our matrix elements can be used at all 
e+e" colliders. 

As a byproduct of the main calculation we have also examined the dominant con- 
tribution, e^e~ QQggi, to the production of a single hard photon in association with 
four jets. Both initial state and final state radiation have been included at tree level, 
without ISR resummation. 

We have computed all matrix elements both in the formalism of [|], ^ and in that of 
||T0| using the former for the actual calculation. The relevant formulae are summarized 
in the following section. Both methods can be easily implemented in a small set of 
nested subroutines. This however results in a computer program which is far too slow 
for the present calculation. We have instead used the symbolic package Mathematica 
| rT| to write down the Fortran expression for each helicity amplitude. The individual 
Z-functions are computed first and then combined in larger structures that describe 
chunks of increasing size of the Feynman diagrams shown in fig. 1-2. Every sub-diagram 
is saved and then used several times. With this procedure we have produced a rather 
large piece of code, which however runs quite fast, and therefeore can be used in high 
statistics Montecarlo runs. As an example the program for qqggg production is about 
24,000 lines long, but requires only about 5 x 10"^ seconds to evaluate on a Vaxstation 
4000/90. 

The amplitudes have been checked for gauge invariance. We have used Mz = 91.1 
GeV, Tz = 2.5 GeV, sm^iOw) = -23, mt = 5. GeV, a^m = 1/128 and a, = .115 in the 
numerical part of our work. 



Calculation 

We present the helicity amplitudes in the formahsm of [|], The two basic functions 
which are needed in writing the amphtudes are: 

Z{pi, Xi- pj, Xj]pk, Xk, Pi, A/; cr, cl; c^, c^) = (1) 
[u{pi, Xi)T''u{pj, Xj)][u{pk, Xk)T'^u{pi, Xi)] 

where 

T^'^'^ = Y{c^r^Pr + cI^Pl) (2) 

and 

X{pi, Xi] Pj;pk, Xk) = [u{pi, Xi)^j u{pk, Xk)] (3) 

Their expressions in terms of spinor products s{i,j) = u{pi, +)u{pj, — ) = t*{j, i), r] and 
/i functions can be found in Table I and II. 

A convenient, though non unique, choice for the spinor products and for the func- 
tions r] and /i is: 

= ipf + 'Pt) (2^1)"^ - (P? + (}^)"' (4) 
r,(i) = (2(p»-pr))'" ''W = =^^ 

where the sign +(— ) refers to a particle (antiparticle) of mass mj. 

The polarization vectors of a gluon of momentum pi can be written as 

e'^iPi, X) = Ni[u{p,, X)^^u{pa(,v), A)] (6) 

where Pa{i) is an auxiliary massless momentum not parallel to Pi and the normalization 
factor is Ni = [4(pa(i) ■ Pi)]~^^'^- AH amplitudes are independent of the Pa(i)'s. 
There are fifty Feynman diagrams contributing to 

e~{pi, A) + e+(p2, -A) q{pG, Xq) + q{p7, Xj) + g{p3, A3) + g{p4, A4) + g{p5, A5) (7) 

Some of them are shown in fig.l. All others can be obtained through permutations of 
the gluon labels. The two diagrams with four-gluon vertices are not shown. Their con- 
tribution is included in the expressions given below for the diagrams with two connected 
three-gluon vertices. 

The amplitude squared for this process is: 

2 1 1 



M 



|A| 7r,7r' 



where vr is a permutation of the indexes (3,4,5), {A} indicates the helicities of all 
external particles, M^^ is the coefficient of the color matrix 



in the full amplitude and C^^/ is the appropriate color factor, C^tt' = Tr{T^(T'^'y). 

The color factors arc given in Table III. To fix our notation, diagrams 1-4 have 
color structure 345, diagrams 5-7 have color structure [3,4]5 where [3,4] is the usual 
commutator, [3,4]=34— 43, diagrams 8-10 have color structure 5[3,4] and diagrams 11- 
12 have color structure [[3, 4], 5]. 

The production of four quarks and one gluon 

e"(Pi, A) + e+(p2, -A) g{p3, A3) + g(p4, A4) + g(p5, A5) + g(p6, Ae) + q{p7, A7) (10) 

is described in lowest order by twenty-four Feynman diagrams if the two quarks have 
different fiavour (case 1) and by forty-eight if their fiavour is the same (case II). Twelve 
of the diagrams for case I are shown in fig.2, all others can be obtained with the simul- 
taneous substitutions q{p5, X5) q{p7,M) and q'(p4, A4) ^ q{pQ,\Q). The additional 
ones needed for identical flavour can be obtained exchanging q{p5, A5) with q{p7, A7) 
including their color quantum numbers and changing the overall sign. 
The amplitude squared for these processes is: 
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in case I and case II respectively. Cij and C'^j are the appropriate colour factors which 
can be extracted from Table IV and V. 

We define the following colour structures: 



n]2 



and 

The color factor is given by: 



(14) 
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in case I and by analogous formulae in case II, possibly with A„,rn B„ 

We adopt the following shorthand notation, where in the left hand side the index 
i stands both for the momentum pi and for the corresponding polarization Aj, and [i] 



stands for the set of four indices associated with a gluon momentum and to his auxihary 
momentum: 

X{i;j;k) = X{pi, Xi;pj;pk, Xk) (17) 

= X{pi,Xi;pj;Pa(i)Ai) (18) 

Zii,j;k,l) = Z{pi,Xi;pj,Xj;pk,Xk;pi,XrAA'AA) (19) 

Zi\i];j,k) = Z{pi, Xi;pa(i), Xi]Pj, ^j]Pk, XkAA'AA) (20) 

[j]) = Z{pi,Xi;Pa(i),Xi;pj,Xj;Pa(j),Xj;l,l;l,l) (21) 

ZeiiJ) = ^(p»,Ai;pj,Aj;p2,A;pi,A;l, 1;1,1) X (22) 
+ Z{i, Xi] j, Xj] 2, A; 1, A; cr, cl, c^r, Cl) x dz^ 

where c?^, dzo are the photon and Z^ propagator functions. The attachment of two 
external gluons to a quark hne through a three-gluon vertex is conveniently described 
with the help of: 

Z2{\i,3];k,l) = Z{[^-[j])x[x{k-3-l)-X{k-i-l)^ (23) 
+ 2(x{[j\-i) X Z{[^-k,l)-X{\i\-j) X Z{\3]-k,l) 



The function that describes the insertion of two connected three-gluon vertices, and 
of the term with the same color structure in the associated four-gluon vertex, on a 
fermion line is: 

^3([«, j], [k]]hrn) = 

4Z([2];/,m) xX([j];z) x (x ([k]; j) + X {[k]- 

+ 2Z{[k]-l,m) X ^Z{\i]-\3])xp.i-pk 

+X([2];j) X (x([j];z) + 2X([j];fc))) (24) 

+ {x{l; k] m) - X{1] i; m) - X{1; j; m)^ x 

Z([^];[j])xX([A:];j)-2^([A:];b1)xX([^];j) 



+ 2X(/;j;m) x Z{\i]-\j]) x [X{[k]- j) + X{[k]-t] 
-2Z{[i\-l,m)xZ{[k]-\j])xp.,-p^ 

The denominator associated with this function is {pi +PjY{pi + Pj +PkY- 

In order to simplify our formulae we make the convention that all repeated indices 
are summed over, and that the sum extends to all possible values, A = ±1, of all internal 



helicities. The last line in eq. 26-49 gives the set of values taken independently by each 



momentum index. Defining the sign factors e{i) as: 



the spinor functions necessary for qqggg production, which correspond to the diagrams 
shown in fig.l are: 

Ml = -Z{[3];6,i)xZ{[4];t,j)xZ{[5];j,k)xZ,{k,7)x6{k) (26) 

Pi = {P6,P3}, Pj = {P&,P3,P4}, Pk = {P1,P2,P7} 

M2 = -Z{[3];6,t)xZ{[A];t,j)xZ,{j,k)xZ{[5y,k,7) (27) 

Pi = {P6,P3}, Pj = {P6,P3,Pa}, Pk = {P5,P7} 

M3 = Z{[3];6,t)xZ,{t,j)xZ{[A];j,k)xZ{[5];k,7) (28) 

Pi = {P6,P3}, Pj = {P4,P5,P7}, Pk = {P5,P7} 

M4 = Zei6,i)xZi[3];i,j)xZ{[A\;j,k)xZ{[5];k,7)xe{i) (29) 

Pi = {P6,Pl,P2}, Pj = {P4,P5,P7}, Pk = {P5,P7} 

Ms - -Z2([3,4];6,2)xZ([5];^,j)xZe(j,7)x£(j) (30) 

Pi = {P6,P3,Pa}, Pj = {Pl,P2,P7} 

Me = -Z,{[3,A];Q,t)y<Ze{t,j)xZ{[5];j,7) (31) 

Pi = {P6,P3,P4}, Pj = {P5,P7} 

Mr = -Z,{Q,t)xZ2{[3,4];^,j)y<Z{[5];j,7)xe{i) (32) 

Pi = {Pl,P2,Pe}, Pj = {P5,P7} 

Ms = -Z{[5]■,6,^)xZ,{[3A]■,^,J)y<Ze{J,7)x6{J) (33) 

Pi = {P6,P5}, Pj = {PUP2,P7} 

Mg = -Z{[5];6,t)xZ,{t,j)xZ2{[3,4];jJ) (34) 
Pi = {Pe,P5}, Pj = {P3,Pa,P7} 

Mio = -Z,{6,t)xZ{[5];t,j)xZ2i[3A];j,r)xs{{) (35) 

Pi = {Pl,P2,P6}, Pj = {P3,P4,P7} 

Mn = -Z3([3,4],[5];6,z)xZe(z,7)x£(z) (36) 

Pi = {Pl,P2,P7} 



M12 = Z,{6,i)xZs{[3,4],[%h7)xe{i) 

Pi = {Pl,P2,P6} 



(37) 



The spinor functions for the diagrams in fig. 2, which describe qqqqg production are: 
Mi3 = Z([3];i,5)xZ(4,i;6,j)xZe(i,7)x£(j) (38) 

Pi^{P3,P5}, Pj ^ {Pl,P2,P7} 

Mi4 = -Z{[3];A,t)xZ{i,5;6,j)xZ,{j,7)xe{j) (39) 

Pi = {P3,P4}, Pj = {Pl,P2,P7} 

Mi5 = -Z([3];i,5)xZ(A,i;j,7)xZ,(6,j)xe(j) (40) 

Pi^{P3,P5}, Pj ^ {Pl,P2,P6} 

Mi6 = Z{[3];4,t)xZ{t,5;j,7)xZ,{6,j)xe{j) (41) 

Pi = {P3,Pa}, Pj = {Pl,P2,P(i} 

M^r = -Z(6,rA,5)xZ{[3];i,j)xZ,{j,7)xe{j) (42) 

Pi = {P4,P5,P6}, Pj = {Pl,P2,P7} 

Mi8 = -Z(6,z;4,5)xZe(z,j)xZ([3];j,7) (43) 

Pi = {P4,P5,P6}, Pi = {P3,P7} 

Mi9 = -Z,(6,z)xZ(z,j;4,5)xZ([3];j,7)x£(z) (44) 

Pi = {Pl,P2,P6}, Pj = {P3,P7} 

M20 = -Z([3];6,i)xZ(i,i;4,5)xZe(i,7)x£(j) (45) 

Pi = {P3,P6}, Pj ^ {Pl,P2,P7} 

M21 = -Z{[3];6,t)xZ,{t,j)xZ{j,7;4,5) (46) 

Pi = {P3,P6}, Pj = {P4,P5,P7} 

M22 = -Z,(6,i)xZ{[3];i,j)xZ(j,7-A,5)xe(i) (47) 

Pi = {Pl,P2,P6}, Pj = {P4,P5,P7} 



M23 = Ze(6,^)x (^2(^X([3];4) + X([3];5)j xZ(z,7;4,5) 

-2X(4;3;5)xZ([3];z,7) (48) 

- (X(i; 4; 7) + X(i; 5; 7) - X{i; 3; 7)) x Z{[3] ; 4, 5)) x e{i) 

Pi = {Pl,P2,P6} 

M24 = -Ze(i,7) X (2(x([3];4)+X([3];5)) xZ(6,i;4,5) 

-2X(4;3;5) X Z([3];6,i) (49) 

- (^X{6; 4; i) + X{6; 5; i) - X(6; 3; i)^ x Z([3]; 4, 5)) x e{i) 

Pi = {Pl,P2,P7} 

From the examples above the one can to construct the full amplitudes for 5j and 4j7 
production with little effort. 



Results 



Our results are presented in fig.3 and 4. In fig. 3 we show the total cross section for 
five-jet production as a function of Ucut in the JADE scheme [|12| and in the DURHAM 



scheme [0]. We have neglected the mass of the quarks of the first two generations, and 



have summed over all their contributions. The cross section for four 6-quark plus one 
gluon-jet production is given by the dotted lines. The dashed lines represent the cross 
section for the production of two 6-quarks plus any combination of three light-parton 
jets. Finally the sum of all contributions to five-jet production including only light, 
strongly interacting, particles is shown by the continuous line. It is to be noticed that, 
at least in the JADE scheme in which neglecting the mass of the light hadrons which 
are effectively detected, represents an actual invariant mass, there is a natural lower 
limit to the allowed range in Ucut for 6-quark production, which corresponds to the 
invariant mass of the 6, ~ (mb/M^o)^ ^ 3 x 10~^. At lower ycut one begins to 
resolve the decay products of the heavy quark. 

In fig.3 we also give, in the JADE scheme, the cross section for e+e~ Q^QQl-, 
summing over five fiavours. In this case we have taken for the c-quark rrtc = 1.7 GeV. 
Additional cuts are necessary in order to screen our results from the collinear singular- 
ities due to initial state radiation. We have required p] > 5 GeV and | cos 6*^ |< .72; 
these cuts are similar to the selection criteria adopted by the LEP collaborations. 

When our results for five-jet production are compared with the data presented by 
ALEPH III] and OPAL it is clear that the absolute normalization is about a factor 
of five too small. The simplest explanation for this discrepancy is our choice for a^. 
In fact we have used = .115 which corresponds to = M^o with A;^ = 200 MeV 
with five active fiavours in the standard formula: 



«s(Q) 



&oiog(gvA2 



feiiog(iog(gVA^)) 

6glog(gVA2) 



(50) 



The analysis of shape variables and jet rates to 0{al) has shown that, in order to get 
agreement between the data and the theoretical predictions, the scale of the strong 
coupling constant has to be chosen to be Q = x^M^o, with ~ 0.1 0. It has later 
been shown that when the relevant logarithms are properly resummed fl^ agreement 



is obtained for much larger values of the scale, ^ 1. ||T5[|. It therefore not surprising 
that our tree level expressions require a relatively small scale in order to describe the 
data. A related issue, for processes with massive quarks, is the choice of the scale of 
the running mass. Within our tree-level approach we have i gnored this problem, using 
rrib = 5. GeV. 

The cross section for 263j is quite large, of the order of several hundreds pb at 
Ucut = -001 in the DURHAM scheme and ycut = -005 in the JADE scheme, and it 
remains sizable out to ycut ~ -01 and to ycut ~ -02 respectively. Over this limited 
range (j{ycut) decreases by more than one order of magnitude. This behaviour should 
be clearly observable with an efficiency of order .3 for detecting each 6-quark. It will 
be more difficult to detect the production of 46j, whose cross section is only a few ph 
at the lowest values of ycut considered in this paper. 

In fig. 4 we present the cross section ratios a{2b3g) / a{2d3g) (continuous line), a{2u2hg) / a{2u2dg) 
(dashed line) and a{2d2hg) / a{2d2sg) (dotted line) in the two recombination schemes 



as a function of ycut- These curves confirm our previous conclusions that mass effects 
increase with the number of final state light partons. The ratio for the dominant 2q3g 
production process is equal to .58 at ycut = -001 in the DURHAM scheme and to .67 at 
Vcut = -005 in the JADE scheme. It is even smaller for the processes with four quark 
jets in the final state. This corresponds to a 6-=-8% decrease in the predictions for the 
total five-jet cross section. 

An estimate of the effects due to the charm mass can be obtained from the ratio 
(T{2c2g'~f)/a{2u2g^f) which is equal to .93 at ycut = -005 in the JADE scheme. The 
corresponding value for a{2b2gj)/a{2d2g'j) is .63, slightly smaller than in the case of 
2q3g production. 

Five-jet production will be visible also at higher energies. At ^/s — 200 GeV we 

obtain (T(e+e~ — >■ dd3g) = 9.8 x 10~^ pb and a{e~^e~ bb3g) = 8.8 x 10~^ pb, showing 
that the mass of the b decreases the cross section by 10% even at the highest energy 
envisaged for LEP II. 

Conclusions 

The cross section for the production of 6-quarks in five-jet events has been computed 
using helicity amplitudes and the Z-function formalism. The differences with the 
corresponding massless rates and the total cross sections for 2b3j and ibj production 
have been studied at the Z° peak. Mass effects result in a 30-=-40% reduction on 
individual cross sections and in a 6-i-8% reduction of the total five jet production. The 
production of a hard photon in association with four jets has been investigated. 
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Table Captions 



table I The Z-functions for all independent helicity combinations in terms of the func- 
tions SjtjTj and II defined in the text. The remaining Z-functions can be obatined 
changing sign to the helicities and exchanging + with — , s with t and R with L. 

table II The X-functions for the two independent helicity combinations in terms of 
the functions s, t, rj and /i defined in the text. The remaining X-functions can be 
obatined changing sign to the helicities and exchanging + with — , s with t and 
R with L. 

table III The color factors, multiphed by 9, for e+e~ — > qqggg. 

table IV The color factors, multiplied by 9, for e^e~ ?i9i5'2^2- 

table V The additional color factors, multiplied by 9, needed for quarks of identical 
flavour e"'"e~ —>■ qqqq. 



Figure Captions 



Fig. 1 Representative diagrams contributing to e^e~ qqggg- All others diagrams 
can be obtained through permutations of the gluon labels. The contribution of 
the two diagrams involving a four-gluon vertex, which are not shown, are inluded 
in the analytical expression given in the text for the diagrams with two connected 
three-gluon vertices. 

Fig. 2 Representative diagrams contributing to e'^e~ QQQQ- The remaining dia- 
grams can be obtained through simultaneous substitutions q{p5, A5) ^ QiPr, A7) 
and q{p4, X4) q{pe,\e). The additional ones needed for identical flavour can 
be obtained exchanging q{p5, A5) with q{p7, A7) including their color quantum 
numbers and changing the overall sign. 

Fig. 3 Total cross section as a function of i/cut in the JADE and in the DURHAM 
scheme for four 6-quarks plus one gluon-jet production (dotted line), for the pro- 
duction of two 6-quarks plus any combination of three light-parton jets (dashed 
line), and for the sum of all contributions to fiv-jet production including only 
light partons (continuous line). 

Fig. 4 Ratio of massive to massless cross sections for a{2b3g)/a{2d3g) (continuous 
line), a{2u2hg) / a(2u2dg) (dashed line) and a{2d2bg) / (x{2d2sg) (dotted line) in 
the JADE and DURHAM recombination schemes as a function of i/cut- 
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